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'nI" ■ Abstract 

^P I Exact solutions of the Caldeira-Leggett Master equation for the reduced density 

matrix {x'\p{t)\x) for a free particle and for a harmonic oscillator system coupled to 
Q^ ■ a heat bath of oscillators are obtained for arbitrary initial conditions. The solutions 

^^ . prove that the Fourier transform of p{t) with respect to ^^ ^ ' factorizes exactly into 

i-G I a part depending linearly on p[0) and a part independent of it. The theorem yields 

the exact initial state dependence of p{t) and its eventual diagonalization in the energy 
basis. 
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I. Introduction — Feynman and Vernon|l|] pioneered the 'influence functional technique' 
in which a quantum system of interest and its environment are represented by a single 
Hamiltonian and the time development of the reduced density matrix p of the system is 
computed by tracing out the degrees of freedom of the environment. Using this technique 
on the total Hamiltonian 
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in which the environment is modeled as a thermal bath of oscillators with co-ordinate co- 
ordinate coupling to the system coordinate x, Caldeira and Leggett0 derived the following 
high-temperature master equation for the reduced density operator: 
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where Hji is the renormalized system Hamiltonian, 7 is the relaxation rate and 

D = Sm'-fkBT, 
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where ks is the Boltzmann constant and T is the temperature of the bath. The same master 
equation has also been derived from completely different approaches and approximations p, 
3. In particular, for 

2m 
Agarwal's approximations yield Eq(2) with 
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which has been used even at zero temperature by some authors^]. In that case, at all 
temperatures, 

D > Am-fhu. (7) 

For an open system, coupling with this environment may lead to near diagonalization of the 
reduced density matrix in some preferred basis. Such decoherence has obvious conceptual 
interest in understanding the transition from quantum to classical behaviour 0. In addition, 
advances of technology have increased the propsects of tests of decoherence by producing 
superpositions of macroscopically distinguishable states]^. On the other hand, maintenance 
of quantum coherence is crucial to the success of quantum computation, cryptography and 
teleportation^]. 

The purpose of the present work is a rigorous study of quantum decoherence by means 
of exact solutions of the Caldeira-Leggett equation. Some interesting questions which have 
been stimulated by the work of Zurek, Habib and PazH in the weak coupling approximation 



are: which initial states are least susceptible to loss of quantum coherence, and which are the 
likely end states. We are able to give precise answers to these questions without making the 
weak coupling approximation due to a factorization property of the initial state dependence 
of the exact solutions of the master equation. We begin by proving the factorization theorem. 



II. Exact Solution of the Master Equation 

A. Oscillator case — In the position representation we denote 

{x'\p{t)\x) = p{R,r,t), 

where 
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The master equation then becomes: 
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A Fourier transform with respect to R reduces it to a first order partial differential equation. 
Defining 
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To integrate this by Lagrange's method (of characteristics) note that on a curve K 
K{s),t = t{s),r = r{s), the equation becomes 
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We readily obtain the three integrals 
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f/± = C±, Us = Cs, 
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where C± and C3 are integration constants and 
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and A± are constants defined by 






If {K,r,t) and (i^', r', 0) are points on tlie curve (15), Eqs. (16) yield K',r' in terms of U± 

Eq. (17) tiien yields the general solution, 

p{K,r,t)=p{K',r',0)exp{aZ), (20) 
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Eq. (20) is the factorization theorem : the exact solution at an arbitray time is equal to the 
initial reduced density operator with shifted arguments {K',r') times the function exp (aZ) 
which is independent of the initial conditions. It may be verified by direct substitution 
that the expression (20) solves Eq.(12). With p{K,r,t) known, a Fourier transform yields 
p{R, r, t) explicitly. 

B. Free Particle Case — Starting from Eq. (10) with cj = 0, or by taking the cj — > limit 

of Eqs(20)-(22), we obtain. 
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p(K,r,t)=p{K',r',0)lexp——-—KH+^— L({r + r')^ + K] (23) 
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where, 

K' = Ky = ^ + {r- ^)e-''^'K (24) 
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The general solution (23) valid for arbitrary initial conditions agrees with the particular 
solution for an initial Gaussian wavepacket obtained earlier in Ref [11] except that 7 should 
be replaced by 47 in the solution quoted there. 

III. Energy Basis Via Decoherence 

A. Oscillator case — The decoherence mechanism during the system-bath interaction is 
known to suppress the off-diagonal elements of the reduced density matrix of the system in 
an appropriate basis, making all information on the system classically interpretable in that 
basis. For simplified models where the self-Hamiltonian of the system has either been ignored 
or considered co-diagonal with the interaction Hamiltonian, the preferred basis in which the 
density matrix becomes nearly diagonal has been believed to be the one that commutes with 
the interaction Hamiltonian. In the Caldeira-Leggett Model studied here, the coupling to 
the bath is a coordinate-coordinate coupling. Consider the solutions, Eq(20) and Eq(23) 
for the oscillator and the free particle, respectively, at long times {'-ft » 1) [12]. For the 
oscillator the solution at t ^ 00 for the overdamped case (7 >> cj) is: 

P(^'-'^)=P(0'0'0)^-P[-T^;^(^' + ^^)]- (25) 

This limiting density matrix is actually independent of the initial density matrix p(0) pro- 
vided that p(0) is normalized, because then 

p(0.0.0) = li^ = l/v^. 
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The final density matrix is not diagonal in the position basis or momentum basis. For 
example, in the position basis we find, 

(a;'|p(oo)|a;) = Aexp(— a;+(a;^ -|- x'^) — 2q;_xx'), (26) 

where 

Since p(oo) is Hermitian, we seek an orthonormal complete set of states |0„) s-t. 

p(oo)|0„) = A„|0„), (28) 

and 

p(oo) = ^A,|0„)(0„|. (29) 
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Using the relation [13], 
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we find the solutions 0„ of (28) to be, 
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wliere Hn{x) are the Hermite Polynomials, and An are normalization constants. We also 
find the eigenvalues 

Srri'yhw Z D — Am'yhw\"' 
D + Arwyhw \D + Arwyhw J 

Thus the final density matrix is diagonal in the basis of energy eigenstates of the oscillator 
(Eq. (31)). The eigenvalues A„ form a convergent geometric series due to the inequality (7) 
and yield Trp(oo) = 1. 

The Wigner function corresponding to the final density matrix is easily seen to be of 
Gaussian form (in agreement with a theorem of Tegmark and Shapiro [14]), 

w{x,p,t = oo) = /_^^(^-||pMI^ + |K''^' 

It yields the position and momentum uncertainties 
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and 

AxAp = > -n. (35) 
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The linear entropy corresponding to the final density matrix is, 

S = Tr(p(oo) - p\oo)) = 1 - ^Hl^ > 0. (36) 

Since the density matrix for t — ^ oo is independent of the initial state, the final uncertainty 
product AxAp and the entropy production (36) are also independent of the initial state. 
At intermediate times the exact density matrix (20) shows a factorized dependence on the 
initial density matrix which we hope to study in detail later. 

B. Free Particle Case — Eq. (23) shows that due to the factor exp(— DA;^t/(16m^7^)) on 
the right-hand side, the density matrix is driven to a diagonal matrix in momentum space 
for t -^ oo. Further, the diagonal elements {K = p' — p = 0) become 

Since momentum being diagonal implies energy being diagonal for a free particle, the energy 
basis emerges in this example too as the preferred basis for t ^ oo. 
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Conclusion — We have obtained exact solutions of the Caldeira-Leggett Master equation 
for the reduced density matrix for an oscillator and for a free particle for arbitrary initial 
conditions in a compact factorizable form. The solutions for both cases studied show that the 
density matrix eventually diagonalizes in the energy basis at long times though the coupling 
to the bath is via position. Our conclusion is in tune with the recent result of Paz and 



Zurek |10| where they show that eigenstates of energy emerge as pointer states, but we do 
not use any weak coupling approximation. The t -^ oo form of our Wigner function agrees 
with a theorem of Tegmark and Shapiro [14]. We hope to study later the intermediate time 
behavior of the density matrix using the factorised exact solution obtained here. 
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